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ABSTRACT 

An  inter'  al  function  Y  assigns  an  interval  Y(x)  *=  I^(x),y(x)J  in  the  ex¬ 
tended  real  number  system  to  each  x  in  its  interval  X  *  [a,b]  of  definition. 

b  b 

The  integral  of  Y  over  ia,b]  is  taken  to  be  the  interval  / Y(x)dx  *  [Jj£(x)dx  , 

a  a 

_b  b 

/y(x)dx],  where  J%(x)dx  is  the  lower  Darboux  integral  of  the  lower  endpoint 
;  a  ~a  , 

—  b 

function  and  /  y(x)dx  is  the  upper  Darboux  integral  of  the  upper  endpoint 

a 

function  y.  As  these  Darboux  integrals  always  exist  in  the  extended  real  number 
system,  it  follows  that  all  interval  functions  are  integrable,  no  matter  how  nasty 
the  endpoint  functions  y_,  y  are.  The  interval  integral  defined  in  this  way  in- 

■  eludes  the  interval  integral  of  R.  E.  Moore  as  the  special  case  that  y  are 
continuous,  and  hence  Riemann  integrable. 

In  addition  to  a  construction  of  the  interval  integral  in  a  form  suitable 
,  for  numerical  approximation,  some  of  its  basic  properties  and  other  implications 
and  applications  of  its  definition  are  presented.  The  theory  of  interval  inte¬ 
gration  given  here  supplies  a  previously  lacking  mathematical  foundation  for  the 

i 

■  numerical  solution  of  integral  equations  by  interval  methods. 

■i 
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SIGNIFICANCE  AND  EXPLANATION 


Solving  physical  problems  with  the  aid  of  digital  computers  often  involves 
performing  hundreds  or  thousands  of  inexact  arithmetic  operations  on  inaccurate 
data.  A  study  of  the  effect  of  these  errors  on  the  final  result  to  determine  if 
it  is  useful  (or  even  meaningful)  is  often  required,  and  may  be  tedious.  By  the 
methods  of  interval  analysis,  the  computation  may  be  carried  out  automatically  on 
intervals  which  contain  the  data  to  yield  intervals  in  which  any  exactly  computed 
result  would  lie.  It  can  also  be  argued  that  the  intervals  calculated  are  actual¬ 
ly  the  solutions  to  certain  problems;  for  example,  the  intervals  of  deflection  of 
a  structure  corresponding  to  intervals  of  loading  forces  may  give  a  better  pic¬ 
ture  of  its  stability  and  safety  than  knowing  its  exact  configuration  under  a 
specific  selection  of  forces. 

Whether  one  interprets  the  obtained  interval  to  contain  the  answer  or  be  the 
answer,  it  is  important  to  be  able  to  compute  interval  solutions  efficiently,  with 
as  little  extraneous  expansion  as  possible.  However,  the  subject  of  interval  anal¬ 
ysis  is  still  in  the  early  stages  of  its  development.  Algebraic  operations  on,  in¬ 
tervals  are  well  understood,  and  can  be  expressed  in  terms  of  the  endpoints  of  the 
intervals  involved.  \In  the  case  of  transformations  of  interval  functions  Y  (Y (x) 
is  an  interval  for  each  x) ,  this  simple  situation  may  not  hold,  except  in  special 
cases.  In  particular ,  the  lack  of  a  theory  of  integration  of  interval  functions 
has  hampered  the  development  of  interval  methods  for  the  solution  of  integral 
equations,  which  constitute  an  important  class  of  mathematical  models  of  real 
.situations.  This  report  provides  such  a  theory  by  defining  integrals  of  interval 
functions  to  be  an  interval  with  the  limits  of  integrals  of  step-functions  respec¬ 
tively  less  and  greater  than  the  integrand  as  endpoints.  As  these  lower  and  up¬ 
per  Darboux  integrals  of  basic  integration  theory  always  exist,  it  follows  that 
all  interval  functions  are  integrable.  The  simplicity  of  this  theory  is  due  to 
the  fact  that  intervals,  rather  than  real  numbers  only,  are  accepted  as  values  of 
integrals,  even  if  the  integrand  is  a  real  function  (called  a  degenerate ^interval 
function  in  this  context) Progress  in  mathematics  is  often  made  in  this  way,  by 
introducing  a  new  object  a i  a  solution;  for  example,  imaginary  numbers  as  square 
roots  of  negative  numbers.  -^Some  useful  properties  of  the  interval  integral  de¬ 
fined  in  this  report  are  derived,  and  its  potential  usefulness  in  connection  with 

the  numerical  solution  of  integral  equations  is  indicated.  ...  -  «— 
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summary  lias  with  MRC,  and  not  with  the  authors  of  this  report. 


INTEGRATION  OF  INTERVAL  FUNCTIONS 
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1 •  Intervals  in  the  extended  real  number  system.  In  ordinary  interval 
analysis  [5],  [6] ,  the  term  interval  refers  to  closed  intervals  of  real  numbers, 

(1.1)  X  =  [a,b]  =  {x  j  a  <_  x  <_  b} , 
with  finite  endpoints  a,  b.  The  width 

(1.2)  w (X)  =  w([a,b])  =  b  -  a 


of  an  interval  with  real  endpoints  is  consequently  finite.  To  develop  the  theory 
of  integration  of  interval  functions  given  below,  it  is  convenient  to  use  the  ex¬ 
tended  real  .number  system,  which  includes  the  values  +  «  [3  ] .  Thus,  in  addition 
to  finite  intervals  of  the  form  (1.1)  with  a,  b  finite,  there  will  be  infinite 
intervals  in  the  system  of  one  of  the  following  types: 


(i)  semi-infinite  intervals 

(1.3)  S  =  [a,  +<*>]  ,  Sb  =  [-»,  b]  ,  a,  b  finite; 

cl 

(ii)  the  real  line 

(1.4)  R  =  [-«,  +»] , 

and  (iii)  the  indegenerate  intervals 

(1.5)  S’"  =  [-»,  -®1 ,  S  =  [+”,  +«] . 

(In  what  follows,  "+°°"  will  often  be  written  simply  as  "«>" . ) 


_ All  the  infinjteintervals  will  be  defined  to  be  of  infinite  width,  that  is, 
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(1.6) 


w(S  )  =  w(S  )  ■  w(R)  =  w(S  )  =  w(S^  )  *  +® 
a  +® 

in  the  extended  real  number  system.  This  definition  is  consistent  with  the  type 
of  limiting  process  used  to  define  "improper"  integrals,  that  is, 


(1.7)  sa  " 

Si  oo 

and  hence 


(1.8) 


, \  lim  ,  _  >  lim  ,  . 
w (S  )  =  _  w (S  )  =  (“)  = 

co  a-«°  a  a-*° 


so  that  it  is  reasonable  in  this  sense  to  assign  infinite  widths  to  the  inde- 

—  00 

generate  intervals  S  and  S  . 

oo 

In  what  follows,  a  closed  interval  in  the  extended  real  number  system  will 
be  called  simply  an  interval . 


2 .  Interval  arithmetic.  Interval  arithmetic  ( [  5  ]  »  [  6  1 )  as  defined  for  finite 
interval  may  also  be  performed  in  the  system  of  intervals  in  the  extended  real 
number  system  defined  in  §1  if  suitable  rules  are  adopted  for  computing  with  the 
values  +«>.  In  essence,  these  "rules"  are  a  shorthand  notation  for  the  results 
of  the  types  of  limiting  processes  to  be  encountered  in  the  theory  of  integra¬ 
tion  presented  below.  McShane  ((3  1,  p.  21)  gives  the  following  rules: 

(i)  -<«>  <  a  <  oo  for  every  real  number  a; 

(ii)  “>,a  =  a*00  =  ®  if  0  <  a  <  ®i 

(iii)  °°*a  =  a**  =  -«*  if  a  <  0; 

(iv)  (-<*>)  *a  =  a*  (-»»)  =  -oo  if  0  <  a  <  °°; 


(2.1) 


(v)  (-oo)  »a  =  a*  (-oo)  =oc  if  -oo  ^  a  <  oj 

a  a 


(vi)  - 


(-«) 


=  0  if  a  is  real; 


(vii)  oo  +  a  =  a  +  oo  =  oo  if  a  >  -oo; 
(viii)  -oo  +  a  =  a  +  (-oo)  =  -«  if  a  <  oo. 
(ix)  oo«o  =  0,x  =  (-oo)  »o  =  0*  (-oo)  =  o. 


Thus,  rule  (2.1) (ix)  takes  care  of  the  "indetenninant  form  '0»®,H  which 
can  arise  if  one  of  the  factors  in  a  multiplication  is  am  infinite  interval.  The 
product  of  two  intervals  will  be  defined  to  be 

(2.2)  [a,b]*[c,d]  =  [min{ac,ad,bc,bd} ,max{ac,ad,bc,bd}] 

in  the  extended  real  number  system.  In  ordinary  interval  arithmetic  ([5  ],  p.  9), 

(2.2)  is  a  consequence  of  the  definition  [a,b]*[c,d]  =  {z  |  z  =  x»y,  x€[a,b], 
y€[c,d]}  of  multiplication  of  intervals.  In  the  extended  real  number  system, 
however,  one  has  {z  ]  z  =  x*y,  x€[-l,l]  ,  y€[®,®]}  =  {-®,  0,  “}  by  (2 .1)  (ii)  (iii) 
(ix) ,  and  the  result  is  not  an  interval.  Use  of  the  rule  (2.2)  gives  [-1,1]* 

[®,®]  =  [-oofo>]  ,  which  circumvents  this  problem. 

As  in  ordinary  interval  arithmetic,  division  by  intervals  containing  0  will 
not  be  defined.  The  reciprocal  of  an  interval, 

(2.3)  [c,d]_1  =  (1/d  ,  1/cl,  OjE[c,d]  , 

is  defined  for  all  zero-free  intervals,  with  rule  (2.1) (vi)  used  if  [c,d]  is  an 
infinite  interval.  One  has 

s"1  =  [a,®]"1  =  [0  ,  1/a],  a  >  0;  s”1  -  [®,®]-1  =  [0,0]; 

a  ® 

(Sb)-1  =  [-®,b]_1  =  [1/b  ,  0] ,  b  <  0;  (s'”)_1  =  [-®,-®]"1  -  [0,0]. 

The  indeterminant  form  "•/«>»  will  thus  not  occur  in  the  interval  arithmetic  un¬ 
der  discussion,  as  division  is  defined  by 

(2.5)  [a,b]/[c,d]  *  [a,b]*[c,d]  1,  0jE[c,d]  , 

and  [c,d]  1,  if  it  exists,  will  have  only  finite  endpoints  by  (2.3)  and  (2.4). 

The  indeterminant  form  "®  -  ®"  cam  appear  in  addition  or  subtraction  ac¬ 
cording  to  the  usual  rules  ( [  5  ]  ,  pp.  8-9)  , 


-3- 


(2.6) 


[afb]  +  [c,d]  -  [a  +  c  ,  b  +  d] 


[a,b]  -  (c,d]  *  [a  -  d  ,  b  -  c]  , 

but  only  if  at  least  one  of  the  terms  is  an  indegenerate  interval.  Thus,  an  ad¬ 
ditional  rule  to  augment  the  list  (2.1)  is  needed,  which  is 


(2.7) (x) 


[a,“]  +  =  (a,®]  -  [<*>,«]  =  [-“,*] 

[-°°,b)  +  =  [-®,b]  -  [-<*>,-»]  00  [-00,00]  , 


where  a,b  may  be  finite  or  infinite.  Thus,  rule  (2.7) (x)  assigns  the  value  +°° 
to  oo  —  oo  as  an  upper  endpoint  of  an  interval,  and  -<*>  as  a  lower  endpoint. 

Thus,  the  total  collection  of  rules  for  interval  arithmetic  in  the  system 
of  intervals  defined  over  the  extended  real  numbers  consists  of  (2.1) (i)-(ix) , 
(2.7) (x),  (2.2),  (2.3),  (2.5),  and  (2.6).  The  interval  arithmetic  constructed 
in  this  way  contains  ordinary  interval  arithmetic  on  finite  intervals  ( [  5  ]  ,  [  6  ] ) 
in  the  sense  that  it  gives  the  same  results  for  finite  intervals.  The  operations 
on  infinite  intervals  are  defined  in  such  a  way  as  to  be  convenient  in  the  sequel 
for  the  construction  of  a  theory  of  integration  of  interval  functions.  Other  ex¬ 
tensions  of  interval  arithmetic  to  infinite  intervals  are  possible,  but  will  not 
be  considered  here. 

3.  Interval  functions.  Y  is  said  to  be  an  interval  function  of  x  on  [a,b]  if 
it  assigns  a  nonempty  interval 

(3.1)  Y(x)  =  ty(x)  ,y(x)]  =  {y  |  £(x)  <  y  <  y(x)} 

to  each  x€[a,b] .  The  (extended)  real-valued  functions  £,  y  are  called  the  end¬ 
points  ,  or  boundary  functions  of  Y,  and  the  notation 

(3.2)  Y  =  [£  ,  y] 


will  be  used,  as  well  as  the  alternative  notation 
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(3.3) 


Y(x)  =  ,  y]  (x) 


for  the  interval  (3.1). 

The  interval  function  Y  can  also  be  identified  with  its  graph,  which  is  the 
set  of  points 

(3.4)  Y  =  la,b]*Y(x)  ■  {(x,y)  |x£[a,b],  y€Y(x)} 

in  the  xy-plane.  Geometrically,  the  graph  (3.4)  extends  from  the  "lines"  x  *  a 
on  the  left  to  x  =  b  on  the  right,  and  from  the  "curves"  defined  by  y  *  £(x)  be¬ 
low  to  y  **  y(x)  above  (recall  that  extended  real  values  are  permitted). 

In  the  context  of  interval  functions,  a  real-valued  (or  extended  real  val¬ 
ued)  function  f  is  considered  to  be  the  degenerate  interval  function 

(3.5)  f  =  [f  ,  f]  . 


In  the  extended  real  number  system,  numbers  c  <_  d  exist  such  that  the  graph 
(3.4)  of  Y  is  contained  in  the  rectangle  R  *  [a,b]x[c,d]  *  ((x,y)  |  x€(a,b] ,  y€ 
[c,d]}  in  the  xy-plane,  that  is, 

(3.6)  Y  c  [a,b] x[c,d]  *  R. 


The  set  of  all  rectangles  R  for  which  (3.6)  holds  will  be  denoted  by  R(Y)  or  by 
R,_  ,_,(Y)  if  it  is  desired  to  specify  the  interval  of  definition  [a,b]  of  Y. 

If  [a,b]  is  a  finite  interval,  then  Y  is  said  to  be  finitely  defined.  If 

(3.6)  holds  with  c,d  finite,  then  Y  is  called  a  bounded  interval  function.  A 
bounded  and  finitely  defined  interval  function  is  said  to  be  finite;  the  graph 
of  a  finite  interval  function  is  obviously  contained  in  a  finite  rectangle  R  with 
area  w( [a,b3 ) *w( [c,d] }  ■  (b  -  a)»(d  -  c). 

Definition  3.1.  For 


inf 

x€[a,b] 


{^(x)}. 


d  ■ 


(3.7) 


c 


8up  (7(X)  } 

x€[a,b]lylxU' 


the  interval 


(3.8)  VY[a>b]  =  tc,d] 

is  called  the  vertical  extent  of  the  interval  function  Y  on  [a,bj .  If  the  inter¬ 
val  of  definition  of  Y  is  understood,  then  VY.  .  ,  may  be  abbreviated  as  VY. 

[a,b]  * 

The  rectangle 

(3.9)  R(VY)  =  R ,  (VY)  =  [a,b]xVY. 

[a,bj  [a,b] 

is  the  "smallest"  containing  the  graph  of  Y.  One  has 

(3.10)  R (VY)  =  0  R  ; 

RGr(Y) 

that  is,  R(VY)  is  the  intersection  of  all  rectangles  (3.6)  which  contain  the 
graph  of  Y. 

Vertical  extent  of  an  interval  function,  as  defined  above,  has  the  important 
property  of  being  inclusion  monotone  with  respect  to  the  interval  of  definition 
of  the  interval  function  and  inclusion  of  interval  functions;  Y  c  z  means  that 
the  graph  of  Z  contains  the  graph  of  Y  considered  to  be  point-sets  in  the  xy- 
plane.  More  precisely,  (3.7)  and  the  definition  (3.8)  of  vertical  extent  lead 
directly  to  the  following  result; 

Lemma  3.1.  If  I,  J  are  intervals  on  the  x-axis  with  I  c  j,  then 

(3.11)  VYt  c  VY: 

X  J 

if  Y,  Z  are  interval  functions  on  X  =  [a,b]  such  that  Y  c  z,  then 


4.  Vertical  measure  and  Darboux  sums.  Definition  4.1.  The  interval 


,4-11  Vm,y> 1 

! 

is  called  the  vertical  measure  of  the  interval  function  Y  on  [a,b].  Note  that  . 

this  quantity  is  interval-valued,  and  specifies  the  interval  of  definition  of 
the  interval  function  Y  on  which  its  vertical  extent  VY  is  obtained. 


(The  goal  in  this  paper  is  to  construct  a  theory  of  Riemann-type  integrals 


of  interval  functions.  The  horizontal  measure 

(4.2)  H r-  K1(V)  =  [a,b] *w(7Y  ) 

[a,b]  [a,b] 

of  Y  on  [a,b]  may  be  useful  in  a  Lebesgue-type  integration  theory,  but  this  will 
not  be  pursued  further  here.) 

Remark  4.1.  Vertical  measure  is  inclusion  monotone  with  respect  to  inclu¬ 
sion  of  interval  functions:  If  Y  c  z,  then  Wr  ,  , (Y)  c  w.  ^,(Z).  Horizontal 

[a,b]  [a,b] 

measure  is  inclusion  monotone  with  respect  to  the  interval  of  definition:  If 

I  <=J,  then  H  (Y)  c  H_(Y)  . 

X  J 

Proof:  The  assertions  of  Remark  4.1  follow  immediately  from  Lemma  3.1. 

As  usual,  a  set  of  points  {x  .  x, ,  ...,  x  }  such  that 

0  1  n 

(4.3)  a  *  x„  <  x,  <  . . .  <  x  ,  <  x  =  b 

0  —  1  —  —  n-1  —  n 

defines  a  partition 

(4.4)  A  =  (X.,  X_,  ...,  X  ) 

n  la  n 

of  the  interval  X  =  [a,b] ,  where 


(4.5) 


Xi  —  [x^_^  ,  x^ ],  i  =  1,2, .. . ,n. 


n  n 

X  =  U  X  ,  w(X)  =  l  w(X. ) 
i=l  i=l  1 


Obviously, 

(4.6) 

Definition  4.2.  The  interval 

(4.7) 


I.  Y  =  l  w(X  )-VY  -  l  W  , 

An  lxi_i  • 


(Y) 


is  called  the  Darboux  sum  of  the  interval  function  Y  corresponding  to  the  parti¬ 
tion  A  of  X  *  [a,b] ,  where  VY,  *  VYv  =  VY.  has  been  written  for  brevity, 

n  1  Xi  lxi-l'xiJ 

For 
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in  nondegenerate  intervals  of  the  partition  A^,  and 

(4.16)  s(xi)  =  max{dj  |  £  x^},  i  *  0,l,...,n, 

at  the  partition  points  x. ,  x, ,  . . . ,  x  .  It  follows  that 

01  n 

(4.17)  si(x)  <  ^tx)  <  y(x)  £s^x),  a  £  x  £  b. 

The  properties  of  integrals  of  step-functions  are  well-documented  ((3  ] , 
pp.  54-57),  for  example,  if  s^^  and  are  step- functions  on  an  interval  X,  and 
k  is  a  finite  constant,  then 

(a)  J^s^xjdx  -  k'/^Mdx; 

(b)  /  ( s  (x)  +  s  (x))dx  -  /  s  (x)dx  +  /  s  (x)dx; 

(4.18)  X  1  2  X  1  X  2 

(c)  if  s^fx)  <  s2 (x)  for  all  x  6  X,  then 

j^wd*  £  yxs2(x)dx. 

Furthermore,  if  s(x)  is  a  step-function  on  X,  then  for  each  partition  A  of  X, 

m 

in  1  1 

(4.19)  l  I  s(x)dx  =  /  s(x)dx. 

j-1  Xj  X 

The  integral  of  a  step  function  is  also  invariant  under  translation  ([3  1,  p. 
57). 

The  above  results  may  be  used  to  prove  corresponding  assertions  about  the 
Darboux  sums  (4.7),  taking  into  account  the  differences  between  real  and  inter¬ 
val  arithmetic. 

Theorem  4.1.  If  Y,  Z  are  interval  functions  on  X  =  [a,b]  and  k  is  a  con¬ 


stant,  then 


(4.20) 


(a)  E  k»Y  =  k*E  Y; 

A  A 

n  n 

(b)  Ea  (Y  +  Z)  c  Ea  Y  +  E,  Z; 

A_  A  A 

n  n  n 

(c)  if  Y  c  Z  on  X,  then 


E^  Y  c  E^  Z  (inclusion  monotonicity) . 
n  n 

Proof:  For  finite  k,  (4.20) (a)  follows  directly  from  (4.18) (a);  rule  (2.7) 
(x)  allows  one  to  drop  the  restriction  of  k  to  finite  values.  For  Y  =  ,  y] , 

Z  =  [£  ,  z] ,  the  inequalities 


(4.21) 


,  ,  inf,  •>  inf,  ,  inf,  ■> 

(a)  x  +  z}  >  x  {^}  +  {z} , 

i  i  i 

(b)  <  slpw  +  slp(z). 

i  i  i 


on  the  intervals  X^,  i  =  l,2,...,n,  ((3),  p.  25)  give 

(4.22)  V (Y  +  Z)  c  VY  +  VZv  , 

i  Xi  Xi 


from  which 


(4.23)  Wv  (Y  +  Z)  C  W  (Y)  +  Wv  (Z)  , 

X .  A .  X , 

X  11 

i  =  l,2,...,n,  and  (4.20) (b)  follows.  Finally,  the  inclusion  monotonicity  of 

the  vertical  measure  W  (see  Remark  4.1)  with  respect  to  inclusion  of  interval 
functions  gives  (4.20) (c) .  QED. 

An  analog  of  (4.19)  is  also  available  immediately.  For  m  =  2,  suppose  that 
a  <_  p  <_  b,  and  that 


(4.24) 


is  a  partition  of  X^  = 


(4.25) 


(2)  _ 


(Xll'X12'”,,Xln1) 
(a,p] ,  similarly, 

(X21'X22'”,,X2n2) 
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is  a  partition  of  x2  *  [p,b] .  For  n  *  n1  +  n2 ,  it  follows  that 


(4.26)  An  *  (Xll'X12'**',Xln1'X21","X2n2) 


will  be  a  partition  of  X  ■  [a,b] ,  and 


(4.27) 


z(l)Y  +  E  =  ea  y* 

A  A  n 

nl  n2 


This  may  be  extended  by  induction  to  any  positive  integer  m  >  2. 

A  type  of  mean-value  (or  mean  interval-value)  theorem  holds  for  the  Darboux 
sums  (4.7) . 

Theorem  4.2.  If  X  =  [a,b]  is  a  finite,  nondegenerate  interval,  then  an  in¬ 
terval  Y(A  )  c  VY  exists  for  each  partition  A  of  X  such  that 
n  X  n 


(4.28) 


E  Y  =  w(X) *Y (A  ) . 
An  n 


Proof:  By  (4.7)  and  (4.20) (a). 


(4.29) 


n  w(X. ) 


n  i-1 


say,  where  for  *  w(Xi)/w(X),  i  «  l,2,...,n, 

n  n  n 

(4.30)  r  =  l  a  cL,  s  -  ]>  o^,  l  a  •  1. 

i=l  i*l  i*l 


Thus, 


(4.31) 


min,  ,  _  max,  , 

c  ■  {i>{ci}  i  r  i 


To'diJ  i •  i  “'V  ■ 


and  hence  tr,s]  c  VY  .  Thus,  by  (4.29),  (4.28)  holds  with  Y (A  )  *  [r,s] .  QED. 

x  n 

5.  step- functions  and  Riemann  sums.  For  each  positive  integer  n,  let  denote 
the  set  of  all  step- functions  sn  on  x  having  n  +  1  partition  points  xQ,  x^,  ..., 
disposed  according  to  (4.3) .  Furthermore,  let 


(5.1) 


S  (Y)  »  {s  1  s  €S  ,  s  (x)  <  £(x)  ,  a  <  x  <  b), 
n  — n  - n  n  — n  —  —  — 


S  (Y)  *  {s  |  s  €S  ,  s  (x)  >  y(x),  a  <  x  <  b}. 
n  n  nnn  —  —  — 


The  sets  S  (Y) ,  S  (Y)  are  nonempty,  as  s  e  -»  belongs  to  S  (Y) ,  and  s  =  +®  to 
**  n  — n  n 

S  (Y) .  By  (5.1)  and  (4.18)  (c)  ,  for  each  s  €  S  (Y) ,  s  €  S  (Y) ,  one  has 
n  -n  -n  n  n 


(5.2) 


ib  ,b 

/  s  (x)dx  <_  /  s(x)dx. 


and  consequently 


(5.3) 


•.  *  sSS  </  *. 

— n  — n  a 


inf  rf'- 


{/  s  (x)dx}  <  —  c—  {/  s  (x)dx}  =  d  . 
J  -n  —  s  €S  ‘  n  ~n 

a  n  n  a 


Definition  5.1.  For  each  positive  integer  n,  let  denote  the  set  of  par¬ 
titions  (4.4).  The  interval 


r  y  =  n  e.  y  =  [c  ,  d  ] 

n  A  €0  An  n 

n  n 


is  called  the  Riemann  sum  of  order  n  of  the  interval  function  Y  over  [a,b] . 
Lemma  5.1.  The  interval  £^Y  is  nonempty;  furthermore,  if  m  >  n,  then 


(5.5) 


£  Y  C  £  Y. 
m  n 


Proof:  The  assertion  of  the  nonemptiness  of  the  interval  (5.4)  is  simply 

a  restatement  of  (5.3).  Denoting  the  set  of  Darboux  sums  (4.7)  by  S  ,  if  i  >  n, 

n 


s„  <=  s  , 

n  m 


as  if  £  c  S  ,  then  one  may  take  the  partition  A  defined  by 
An 
n 


(5.7) 


a  *  x0  —  X1  —  *  *  *  —  xn  *  Vl  =  *  Xm  *  b' 


for  which  £  Y  *  E,  Y,  and  thus  £_  Y  €  S  for  m  >  n.  The  inclusion  (5.5)  then 
A  A  Am 

m  n  n 


follows  from  (5.6)  by  the  definition  (5.4).  QED. 
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The  properties  of  Darboux  sums  listed  in  Theorem  4.1  survive  the  intersec¬ 
tion  (5.4)  and  thus  beccme  properties  of  Riemann  sums,  giving  immediately  the 
following  result. 

Theorem  5.1.  If  Y,  Z  are  interval  functions  on  X  =  [a,b]  and  k  is  a  con¬ 
stant,  then 

(a)  E  k»Y  -  k*E  Yi 

n  n 

(b)  I  (Y  +  Z)  c  Z  Y  +  Z  Z; 

(5.8)  n  n  n 

(c)  if  Y  c  z  on  X,  then 

E  Y  e  E  Z  (inclusion  monotonicity) . 
n  n 

The  additivity  of  Riemann  sums  with  respect  to  the  intervals  over  which  they 
are  defined  will  now  be  investigated.  In  order  to  be  definite,  the  notation 


(5.1) 


EY,  ,  .  =  E  Y 
n  [a,b]  n  X 


will  be  used  to  indicate  the  interval  of  summation  X  =  [a,b] .  Suppose  that 
a  <_  p  <_b  and  x^^  =»  [a,p] ,  X2  =  [p,b] .  The  following  results  apply  to  the  expres¬ 
sion  of  the  sum  of  an  interval  function  Y  over  X  in  terms  of  its  sums  over  Xj^ 
and  X^. 

Lemma  5.2.  If  X  *  (a,b]  is  finite  and  nondegenerate,  then 


(5.10) 


’W1' 


Vm"' 


c  W 


[a,b] 


(Y), 


that  is. 


(5.11)  w( (a,p) )  +  w<  [p,bJ)-VY[p  fa)  =w((.,b))"7Y[a  W. 

Proof:  Let  VY[lp)  -  ICj  ,  ^1 ,  -  ICj  .  d,] ,  YY^  -  to  ,  d) . 

Then,  by  definition  of  vertical  extent. 


(5.12) 


c  «  min{c^  ,  Cj) ,  d  ■  max^^  ,  d2> 


For  a  -  w(  [a,p]  )/w(  [a,b] ) ,  one  has  1  a  —  0  311(3  tP»bl  )/w(  [a,b] )  =  1  -  a  >  0. 
Thus, 

(5.13)  w(t«,p]  )-7V[a>  p)  ♦«(tp,b)).VI[pbl  - 

=  w( [a,b])* (ac^  +  (l-a)c2  ,  od^  +  (l-o)d2), 

and,  as 

(5.14)  c  £  ac^  +  (1  -  a)c2  <_  otd^  +  (1  -  a)d2  <_  d, 

by  (5.12),  (5.11)  follows.  QED. 

Theorem  5.2.  If  X  *  [a,b]  is  finite  and  nondegenerate,  then  for  each  pos¬ 
itive  integer  n  >  2, 


(5.15) 


Z  Y,  .  . 

n  (a,b) 


n  {z 


w 


U  Y,  .  +  Z  Yr  c  z  .Yr  .  . 
n  nx  [a,p]  n2  [p,b]J  n-1  [a,b] 


Proof:  The  set  S  of  Darboux  sums  (4.7)  may  be  decomposed  into  two  disjoint 
n 

subsets  for  each  positive  integer  n  >  2:  The  set  ^  of  sums  corresponding  to 

—  n 

partitions  AP  which  have  p  as  a  partition  point,  the  set  of  which  will  be  denoted 
by  ,  and  the  complement  of  relative  to  S  ,  SP  *  S  >>SP,  that  is,  the  set  of 


n  n  n  n 

all  Darboux  sums  corresponding  to  partitions  A  for  which  p  is  not  a  partition 

n 

point.  As  SP  c  S  ,  one  has 
n  n 


(5.16) 


tn!r|«,W  C  n  £A  y[«.bl 
n  n 


By  (4.27) ,  for  Z  Y  €  S  ,  one  can  write 

AP  n 
n 


(5.17) 


ZA  Y.  .  .  =  ZA  Y  ,  +  ZA  Y .  .  . , 

A  [a,b]  A  [a,p]  A„  [p,b] 
n  n,  n„ 


where  n^  +  n2  =  n.  Consequently,  as 
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(5.18) 


n  e4  y. 

AP€V P  An  ta'bl 


n1+n2=n 


{E  Y.  . 
[a»pl 


Z  Y  r 
n2  [p,b] 


the  first  inclusion  of  (5.15)  follows. 

P 

Now,  consider  a  partition  A  of  X  *  [a,bl  for  n  >  2,  and  let  A  denote 

n-1  ~  n 

the  partition  of  X  obtained  by  adding  the  point  p  to  the  set  {xQ,x^, . . . . 
Either  p  =  x.  for  sane  i,  0  £  i  ^  n-1,  in  which  case 

(5.19)  Z .  Y.  .  .  -  Z  Y.  .  .  =  E  Y  .  +  E  Y  ., 

A  ,  [a,b]  p  [a,b]  An  [a,p]  [p,b] 

n  1  2 

n,  +  n„  =  n,  or  a  nondegenerate  interval  X.  =  [x.  ,  ,x.)  c  A  exists  such  that 
12  i  x  n-x 

x.  ,  <  p  <  x. .  As 
1-1  1 


(5.20)  «([*i.1.Pl)*W,Vi,p) 


by  Lemma  5.2,  one  has 
(5.21) 


\pY[a,bJ  C  \_lY[a'bl 


in  this  case.  Thus,  as  {E  Y  .  .}  =  SP, 

Ap  [a,b]  n 

n 


(5.22) 


(1  E  Y.  .  .  c 
Aptf>  fiP  l»pbl  cP.  An-1 

n  n  n  n-1  n-1 


0  E*  Y[a,b]  =  En-lY[a,b] 


by  (5.20)  and  (5.21),  and  the  second  inclusion  in  (5.15)  now  follows  fran  (5.18). 
QED. 

A  mean  interval-value  theorem  also  holds  for  Riemann  sums. 

Theorem  5.3.  if  X  ■  [a,b]  is  finite  and  nondegenerate,  then  for  each  posi¬ 
tive  integer  n,  an  interval  7  c  7Y  exists  such  that 

n  a 

(5.23)  EnY-w(X)*Yn. 
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Proof:  As  before,  let  V  denote  the  set  of  all  partitions  A  for  each  pos- 

n  n 

itive  integer  n.  Then,  by  (4.28) , 


(5.24) 


e  y 

n 


n  e  v 

A  €0  n 
n  n 


w(x) •  n  y(A  ), 

An  &  n 

n  n 


so  that  (5.23)  holds  with 


(5.25)  y  =  n  y(A  )  c  vy  . 

n  A  €D  n  X 

n  n 

as  each  Y(A  )  c  VY.  QED. 
n  X 

6.  Interval  integrals. 

Definition  6.1  (The  Interval  Integral) .  If  Y  is  an  interval  function  de¬ 
fined  on  X  *  [a,b] ,  then  the  interval  integral  of  Y  over  [a,b]  is  defined  to  be 
the  interval 

b  00 

(6.1)  Jy(x)dx  =  Jx Y(x)dx  -  n  EnY. 

a  n=l 

As  usual,  Y  is  said  to  be  integrable  over  X  if  its  interval  integral  (6.1)  is 
defined. 

Remark  6.1.  By  Lemma  5.1,  the  interval  integral  (6.1)  is  a  nonempty  closed 
interval,  as  it  is  the  intersection  of  a  (countable)  collection  of  nested  closed 
intervals . 

Remark  6.2.  An  equivalent  definition  of  the  interval  integral  (6.1)  is 

(6.2)  Jx Y(x)dx  «  [/xZ<x)dx  ,  7xy(x)dx], 

where,  for  the  sets  of  step- functions 

00  00 

(6.3)  S  »  US,  S  *  US, 

—  -i- n  .  n 

n-1  n=l 

the  lower  Darboux  integral  of  ^  over  X  *  [a,b]  is  defined  to  be  ([31,  p.  57)  , 

(6.4)  £xz(x)dx  =  g£s^x®(x)dx^' 

and  similarly,  the  upper  Darboux  integral  of  y  over  X  »  [a,b]  is 
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(6.5) 


/x y(x)dx  -  j^|{/^s(x)dx>, 


The  set  £  defined  in  (6.3)  is  the  set  of  all  step- functions  bounded  above  by 
X,;  similarly,  S  is  the  set  of  all  step- functions  which  are  greater  than  or  equal 
to  y  at  each  point  of  x.  As  these  sets  are  non-empty  (recall  the  step- functions 
^  =  — »  and  s  =  +») ,  the  Darboux  integrals  (6.4)  and  (6.5)  always  exist,  no  matter 
how  nasty  the  functions  y_,  y  are  from  the  standpoint  of  ordinary  integration  the¬ 
ory.  This  observation  furnishes  the  following  result. 

Theorem  6.1  (Theory  of  Interval  Integration).  If  Y  is  an  interval  function 
defined  on  X  *  [a,b] ,  then  its  interval  integral  (6.1)  over  [a,b]  exists. 

In  other  words,  all  interval  functions  are  integrable  (in  the  sense  of  in¬ 
terval  integration) .  The  simplicity  of  this  theory  is  due  to  the  fact  that  inter¬ 
vals  are  accepted  as  values  of  integrals,  including  the  case  that  the  integrand 
is  degenerate  (i.e.,  a  single  real  function).  The  requirement  that  the  integral 
of  a  real  function  be  a  real  number  rather  than  a  possibly  nondegenerate  interval 
leads  to  numerous  difficulties  and  correspondingly  rich  theories  of  integration 
(as  elucidated  in  [  3) ,  for  example),  which  constitute  seme  of  the  most  important 
chapters  of  real  analysis.  By  the  introduction  of  interval  values  for  integrals, 
these  difficulties  are  resolved,  and  the  operation  of  integration  is  extended  to 
all  functions,  interval  or  real.  This  is  analogous  to  the  way  that  the  introduc¬ 
tion  of  complex  numbers  extends  the  operation  of  root  extraction  to  all  numbers, 
complex  or  real.  However,  just  as  complex  analysis  does  not  supercede  real  anal¬ 
ysis,  it  is  to  be  ejected  that  interval  analysis  will  develop  as  a  complementary, 
rather  than  a  competitive,  discipline  to  real  analysis. 

Some  implications  of  the  definitions  of  the  interval  integral  given  above, 
and  some  basic  properties  of  interval  integrals  will  now  be  investigated. 

Remark  6.3.  If  £,  y  are  Riemann  integrable  on  [a,b] ,  then 

b  b  b 

(6.6)  / Y(x)dx  -  [(R)/^(x)dx  ,  (R)/y(x)dx) , 


isy  . . . 


in  terms  of  the  Riemann  integrals  of  the  lower  and  upper  endpoint  functions. 

This  follows  from  (6.2)  and  the  definition  of  a  Riemann  integrable  function 
( t  7  1 »  p.  57)  as  one  with  equal  upper  and  lower  Darboux  integrals;  it  Riemann  in¬ 
tegral  is  taken  to  be  this  common  value,  so  that  if  y  is  a  Riemann  integrable 
function  on  X  *  [a,b] ,  then  its  Riemann  integral  is 

b 

(6.7)  (R)/y(x)dx  =  £  v(x)dx  *  /  y(x)dx. 

a 

Remark  6.4.  In  case  y  are  continuous  on  [a,b] ,  then  the  construction  of 

the  interval  integral  of  Y  may  be  simplified  by  taking  only  the  equidistant  parti- 

tions  A  defined  by  the  points 
- n 

k 

(6.8)  x^  =  a  +  — •  (b  -  a)  ,  k  =  0,l,...,n, 

for  each  positive  integer  n,  so  that  w(x^)  =  and  hence  by  (4.2o)(a), 

i  n 

(6.9)  Ta  Y  =  T  VY.  . 

A  n  .L,  i 
n  i=l 

Here,  the  formation  of  the  Riemann  sums  E  Y  can  be  skipped,  and  the  interval  in- 

n 

tegral  is  given  by 

b  <*>  b  b 

(6.10)  / Y(x)dx  =  0  Z_  Y  =  ((R)/^(x)dx  ,  (R) fy (x)dx] , 

a  n=l  An  a  a 

((31,  pp.  58-59),  as  continuous  functions  are  Riemann  integrable. 

The  interval  integral  (6.10)  is  the  one  proposed  by  R.  E.  Moore  for  contin¬ 
uous  interval  functions  (15  1,  Chapter  8;  [6  1,  pp.  50-56),  as  the  endpoint  func¬ 
tions  of  a  continuous  interval  function  are  necessarily  continuous  (15  ],  p.  18; 

[6  1,  p.  33).  Of  course,  even  in  the  case  y  are  continuous,  one  may  be  able 

to  find  a  partition  A  of  [a,b]  other  than  A  such  that  E.  Y  is  properly  contained 

n  n  A 

n 

in  the  Darboux  sum  (6.9),  and  hence  provides  a  "more  accurate"  approximation  to 
the  interval  integral  than  given  by  use  of  the  equidistant  partition.  Seme  addi¬ 
tional  remarks  about  the  numerical  approximation  of  interval  integrals  will  be 
made  later. 

Some  basic  properties  of  interval  integrals  come  directly  from  the  properties 
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of  the  corrouponding  lUemann  sums  (5.4)  which  hold  under  the  intersections  in 
(6.1).  Thus,  from  Theorem  5.1,  one  has  the  following  result. 

Theorem  6.2.  If  Y,  Z  are  interval  functions  on  X  ■  (a,b)  and  k  is  a  con- 
atant,  then 

b  b 

(a)  /k*Y(x)dx  -  k*/y(x)dx; 

a  a 

b  b  b 

(b)  / (Y (x)  +  Z (x) ) dx  c  / Y(x)dx  +  / Z(x)dxj 

a  a  a 

(6.11) 

(c)  if  Y  c  Z  on  X,  then 

b  b 

/ Y(x)dx  c  / Z(x)dx  (inclusion  monotonicity). 

a  a 

By  taking  intersections  over  all  positive  integers  n  of  the  axprassions  in 
(5.15),  one  gets  immediately i 

Theorem  6.3.  If  Y  is  an  interval  function  defined  on  a  finite,  nondegener- 
ate  interval  X  ■  [a,b] ,  and  p  is  such  that  a  £  p  <_  b,  than 

•  p  b  b 

(6.12)  /y(x)dx  +  /y(x) dx  •  /y(x) dx. 

ape 

Similarly,  Theorem  5.3  furnishes  the  following  mean  interval-value  theorem 
for  interval  integrals. 

Theorem  6.4.  if  Y  is  defined  on  a  finite,  nondegenerate  interval  X  ■  (a,b] , 
then  an  interval  Y  c  7y^  exists  such  that 

b  _ 

(6.13)  /y(x)dx  -  w( [a,b] ) *Y. 

a 

Proof)  Taking  intersections  of  both  sides  of  (5.23)  over  all  positive  in¬ 
tegers  n  gives  (6.13)  with 

00 

(6.14)  Y  -  n  7  C  7Y  . 

n«l 

QED. 

Theorem  6.4  is  useful  in  connection  with  properties  of  indefinite  Integrals. 
Definition  6.2.  The  interval  function 
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(6.15) 


Illx)  =  / Y(t)dt 


is  called  the  indefinite  integral  of  the  interval  function  Y  over  [a  ,  x]  for 
x  ^  a.  (IbY(x)  is  similarly  defined  over  [x  ,  b]  for  x  <_  b.) 

Theorem  6.5.  If  Y  is  a  bounded  interval  function  on  [a,b] ,  then  I  Y (x)  is 

cl 

a  continuous  interval  function  at  any  p  €  [a,b] . 

Proof:  Suppose  that  VY.  =  [c,d] ,  and  take,  as  usual, 

Dj 


(6.16) 


l7Y[a,b]  I  =  max{  H'  ldl>» 


which  is  finite  by  hypothesis.  For  a<p<  b,  a  <  x  <  p,  Y,  ,  c  VY,  ex- 

—  [x,p]  [a,b] 

ists  such  that 


(6.17) 


Xf  *  __ 

/ Y(x)dx  »  /y(t)dt  +  w([x,p])»Y[x 


by  Theorems  6.3  and  6.4;  likewise,  an  interval  Yr  .  cVY.  .  ,  exists  such  that 

[p,x]  [a,b] 


(6.18) 


*  F  _ 

/ Y(t)dt  =  / Y(x)dx  +  w([x,p])*Y[p 


for  p  <  x  <  b.  Given  any  e  >  0,  for  6  =  e/|?Y  . | ,  the  endpoints  of  I  Y(p) 

l&fl)]  3. 

thus  differ  from  the  endpoints  of  I  Y(x)  by  less  than  e  for  |x  -  p|  <  6.  Conti- 

cl 

nuity  of  I  Y(x)  from  the  right  at  x  *  a  and  from  the  left  at  x  *  b  is  obtained 

from  (6.18)  and  (6.17),  respectively,  as  I  Y(a)  =  0  by  Theorem  6.4.  QED. 

£1 

Indefinite  integrals  also  exhibit  a  type  of  differentiability  if  the  limits 


(6.19) 


I’Y(x)  *=  - — -  «T  Y(x)  *  Y 

-  '  ptxw([p,xl)  V  p+xx[p,x] 


(6.20) 


I'Y(x)  «  - — -  »I^Y(x)  *  Y 

+  w  q+x  w ( [x  ,  q] )  q+x  *  [x,q] 


exist  and  are  equal,  where  Y.  and  Y  .  are  the  intervals  defined  in  Theorem 

[PfX]  [x#qj 

6.4,  and  x  lies  interior  to  the  interval  of  definition  [a,b]  of  Y.  (One-sided 
derivatives  at  x  -  a  and  x  *  b  are  defined  by  (6.20)  and  (6.19),  respectively.) 


Definition  6.3.  If  the  limits  I'Y(x)  and  I^Y(x)  exist  and  are  equal,  then 
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(6.21) 


TY  •  I^Y(x)  -  IJY(x) 
io  called  the  derivative  of  the  indefinite  integral  of  Y  at  x. 

The  following  theorem  gives  a  condition  under  which  the  derivative  of  an 
indofinita  interval  integral  is  equal  to  its  integrand. 

Theorem  G.6.  If  Y  is  a  continuous  interval  function  on  (A,b),  then  its  in¬ 
definite  integral  is  differentiable,  and 


I'Y  -  Y (x) ,  a  <  x  <  b, 

(6.22)  x 

I_JY(a)  -  Y (a) ,  TY(b)  ■  Y(b) . 

Proof  i  Let,  for  example,  Y^  ^  ■  (z_(p)  ,  z(p))  for  p  <  x.  Considering  the 
upper  endpoint  function  y  of  Y  as  an  interval  function,  it  follows  from  (4.31), 
(3.25),  and  (6.14)  that  z  €  y,  . .  As  y  is  continuous  if  Y  is  a  continuous  in- 

torval  function, 


6.23) 


imilarly,  ^  z(p)  - 


lim  .  ,  lim  — ,  . 

ptx  ’'IP, xl  ■  •  p+x  lW- 


v(v\ 


_ 4.1.. 4. 


(6.24) 


I'Y(x)  ■  Y(x) ,  a  <  x  <  b. 


In  the  same  way,  one  has 

(6.25)  inr(x)  -  Y (x) ,  a  <  x  <  b, 

which  establishes  (6.22).  QED. 

7.  Relationships  between  interval,  Riemann,  and  Lebesque  integrals  of  real 
functions.  Ordinarily,  no  distinction  will  be  mads  between  a  real  function  y  and 
tho  corresponding  degenerate  interval  function  [y]  ■  (y,y)  having  equal  upper  and 
lower  endpoint  functions.  It  is  convenient,  however,  to  distinguish  between  pos¬ 
sible  integrals  of  y  over  an  interval  X  ■  (a,b] .  The  notation 

b  b  b 

(7.1)  / y(x)dx,  (L)/y(x)dx,  (R)/y(x)dx, 

a  a  a 

will  be  used  to  denote  respectively  the  interval  integral  of  y  as  a  degenerate 
interval  function  (which  integral  always  sxists) ,  the  Lebesgue  integral  of  y  if 
y  is  Lebesgue  integrable  over  [a,b],  and  finally,  the  Riemann  Integral  of^y  if 
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it  OXiatS. 


Remark  7.1.  The  intogral  of  a  dogonoratc  intorval  function  y  i*  »  deganer- 
ato  intorval,  that  is, 


(7.2) 


/y (x)dx  -  [r  ,  r) , 


if  and  only  if  the  real  function  y  ia  Riemann  lntegrable  ovar  [a,b] ,  ao  that 


(7.3) 


r  ■  (R) / y (x)dx. 

a 


This  follows  directly  frcm  Remark  6.3  and  the  definition  (6.7)  of  the  Rie¬ 
mann  integral. 

Thus,  one  ordinarily  expects  an  interval  integration,  even  of  a  single  func- 
tion,  to  result  in  a  nondegenerate  interval.  For  example,  if  xp  i*  the  charac¬ 
teristic  function  of  the  rationale,  that  is, 

(7.4) 


v  (x)  ■  1  for  x  rational, 
P 


X  (x)  ■  0  for  x  irrational, 
P 


then 


(7.5) 


/x_  (x)dx  -  [0,1], 


as  is  well-known.  On  the  other  hand,  seme  nondegenerate  interval  functions  have 
degenerate  interval  integrals.  Consider  the  function  Y  defined  by 


0  <  X  <  y, 


0, 

[0,1], 

(7.6)  Y(x)  -  <  1, 

[1/2], 

.  2, 

i.e.,  Y  is  an  interval  step-function,  which  includes  the  "risers"  as  well  as  the 
"treads".  For  this  function, 


3  <  X  <  f' 1 

x  -  f, 

f  <  x  <  1, 


(7.7) 


/ Y(x)dx  -  [1  ,  1] , 
0 
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au  the  lower  and  uppor  boundary  functions  of  V  have  equal  (Riemann)  integral*. 

Any  interval  function  Y  may  bo  intorprotod,  of  course,  as  a  Bot  of  functions, 
that  is, 

(7<8>  V  =  {y  |yjx)  <  y(x)  <  y(x),  a  <  x  <  b). 

If  Y  is  degenerate,  then  tho  sot  (7.8)  consists  of  only  the  single  function  y  ■ 

Z  “V-  Otherwise,  Ywill  contain  a  number  of  functions,  among  which  there  may 

be  subseta  with  certain  distinguishing  properties  (continuity,  differentiability, 

monotonicity,  etc.).  For  the  discussion  of  integration,  the  following  subsets 

of  functions  will  be  singled  out  for  special  mention. 

Definition  7.1.  If  Y  is  an  interval  function  on  [a,b] ,  then  the  set  of 

Lebesgue  (Riemann)  integrable  functions  y  £  Y  will  bo  called  the  Lebeague  (Ria- 

mann)  core  of  Y,  and  will  be  denoted  bv  C. (Y)  (C_ (Y) ) . 

L  R 

One  has  C  (Y)  c  c  (Y)  always,  but  these  sets  may,  of  course,  be  empty.  For 
R  L 

example,  if  M  is  a  subset  of  [0,1]  which  is  not  measurable  in  the  sense  of  Lebesgue 
then  its  characteristic  function  x»,  is  a  degenerate  interval  function  with  an  emp- 
ty  Lebesgue  (and  hence  Riemann)  core.  The  characteristic  function  xp  of  the  ra- 

1 

tionals  considered  earlier  (see  (7.4))  provides  an  example  of  a  degenerate  inter¬ 
val  function  with  an  empty  Riemann  core,  but  a  nonempty  Lebesgue  core  (the  func¬ 
tion  x  itself ) . 

P 

Definition  7.2.  The  value  v(C_  (Y) )  (v(C_(Y)))  of  the  Lebesgue  (Riemann)  core 

-  L  R 

of  Y  on  [a,b]"is  defined  by 

b 

r  v(Cl(Y})  -  {r  I  r  -  (L) /y(x)dx,  y  £  CL (Y) > , 

(7.9)  ]  b 

1  v (C_ (Y) )  -  {r  |  r  -  (R) /y(x)dx,  y  €  C_(Y)}, 

R  a  R 

respectively,  provided  that  the  indicated  cores  of  Y  are  nonempty. 

Each  set  v(C. (Y) )  and  v(C„(Y)),  when  nonempty,  are  convex,  that  is,  if  one 
L  R  1 

contains  values  r  ,  r2»  with  r1  <_  r2,  then  it  contains  the  entire  interval  [r^, 
r?] .  This  is  because  if  y1  has  integral  r^  and  y2  has  integral  r2,  then  the 
functions  y  ■  y^  +  0(y2  -  y^  are  all  integrable  for  0  <  6  <  1,  and  have 
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■  ^ 

integrals  equal  to  rQ  -  +  0(r2  -  r^ ,  0  <  0  £  1,  which  is  just  another  ex¬ 

pression  for  the  interval  [r^  ,  r,j]  .  As  a  matter  of  fact,  the  theory  of  Lebesgut- 
integration  [  3  ]  leads  to  the  conclusion  that 

(7.10)  I  (y)  =  v (c  (Y) )  , 

L  L 

if  it  exists,  is  a  closed  interval,  which  will  be  called  the  Lebesgue  subinter¬ 
val  of  the  interval  integral  (6.1)  of  Y  over  [a,b]  .  The  set  v(C  (Y)),  on  the 
other  hand,  is  not  necessarily  closed.  This  is  considered  to  be  a  defect  of  Rie- 
mann  integration,  and  led  to  the  construction  of  the  theory  of  Lebesgue  integra¬ 
tion.  However,  as  v(C  (Y))  is  convex,  then  its  closure, 

R 

(7.11)  I  (Y)  =  v (C  ( Y ) ) 

R  P. 

is  a  closed  interval  which,  if  it  exists,  will  be  called  the  Riemann  subinterval 
of  the  interval  integral  of  Y  over  [a,b] . 

The  purpose  of  the  introduction  of  the  intervals  (7.10)  and  (7.11)  is  to 
provide  some  quantitative  information  about  the  Lebesgue  and  Riemann  cores  of  an 
interval  function  Y  which  measures  its  "integrability"  in  a  certain  fashion.  In 
the  metric  topology  for  intervals  ( [ 5  ) ,  [  6  ] )  ,  the  distance  between  intervals 
[a,b]  and  [c,d]  is  defined  to  be 

(7.12)  d( [a,b] ,  [c,d] )  =  max{ |a  —  c |  ,  | b  —  d|}. 

(In  the  extended  real  number  system,  rule  (2.7)  (x)  is  used  to  resolve  any  inde¬ 
terminant  forms  entering  into  (7.12).) 

Definition  7.3.  For 

b 

(7.13)  I (Y)  =  /y(x)dx, 

a 

if  the  Riemann  core  C  (Y)  of  Y  is  nonempty,  then 

R 

(7.14)  p (Y)  =  d(I  (Y)  ,  I(Y)) 

R 

is  called  the  Riemann  gap  of  the  interval  function  Y  on  [a,b] i  similarly,  if 

CT  (Y)  is  nonempty,  then 

Li 

(7.15)  \ (Y)  -  d(I  (Y)  ,  I  (Y) ) 

Lt 
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is  called  the  Lebesgue  gap  of  Y  on  [a,b] . 


Remark  7.2.  One  has 

(7.16)  A  (Y)  <_  p  (Y) 
in  case  both  numbers  are  defined. 

This  follows  from  the  inclusion  C  (Y)  c  c  (Y) .  If  only  one  of  the  numbers 

R  L 

X (Y) ,  p (Y)  is  defined,  it  will  be  X(Y)  by  the  same  token.  For  the  example  (7.4) 

of  the  degenerate  interval  function  x  >  one  has  X(x  )  *  1,  and  p(x  )  is  not  de- 

P  P  P 

fined. 

Theorem  7.1.  If  the  endpoint  functions  £,  y  are  Riemann  integrable  over 
[a,b] ,  then  A(Y)  *  0;  if  A(Y)  =  0,  then  y  are  Lebesgue  integrable,  and 

b  b  b 

(7.17)  /Y(x)dx  *  t(L)/x(x)dx  ,  (L)/y(x)dx). 

a  a  a 

Proofs  By  Remark  6.3,  the  Riemann  integrability  of  y_,  y  means  that  p  (Y)  =0, 
hence  A (Y)  =0  by  (7.16).  Conversely,  if  A(Y)  »  0,  then  the  integral  I (Y)  is 
finite,  and  bounded  sequences  ,  {y^}  c  Y  of  Lebesgue  integrable  functions  may 
be  found  which  converge  to  £  and  y,  respectively.  It  follows  ([3  ],  p.  81)  that 
and  y  are  Lebesgue  integrable  on  [a,b] ,  and,  as 

(7.18)  ^  (Lj/^fxjdx  =  £(x)dx. 


one  has  that 


b 

(7.19)  /.X£(x)dx  *  (L)/x.(x)dx, 

a 

and  similarly  for  y,  whence  (7.17).  QED. 

Remark  7.3.  If  £  and  y  are  Lebesgue  integrable  on  £a,b) ,  then 

b  _  b 

(7.20)  A  ■  max{ (L)/^(x)dx  -  £^(x)dx  ,/xy(x)dx  -  (L)  fy(x)dx) . 

9l  Si 

This  is  true  because  £  is  the  "smallest"  Lebesgue  integrable  function  con¬ 
tained  in  the  interval  function  Y,  and  y  the  "largest"  in  the  sense  that  for  each 
function  y  €  CL(Y) ,  one  has  ^(x)  <_  y(x)  <_  y(x) ,  a  <_  x  <_  b.  Thus, 
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(7.21) 


b  b 

v(Cl(Y))  -  [(L)/^(x)dx  ,  (L)j7(x)dx] , 
a  a 

from  which  (7.20)  follows  by  (7.12). 

8.  Improper  integrals .  In  ordinary  integration  theory,  an  integral 

b 

(8.1)  r  *  j y(x)dx 

a 

is  said  to  be  improper  if  the  interval  of  integration  [a,b]  is  infinite,  or  if 

its  integrand  is  unbounded  on  X  =  [a,b]  in  the  sense  that  given  any  M  >  0,  there 

exists  a  nondegenerate  subinterval  X  of  X  such  that  ] y (x) [  >  M  for  x  €  X  .  Sup- 

posing  that  y  is  unbounded  at  x  =  a,  that  is,  X  «  [a  ,  a] ,  a  ■  a(M)  >  a,  and  in- 

M 

tegrable  (in  the  sense  of  Riemann)  on  [a  ,  b] ,  one  defines  the  improper  Riemann 
integral  of  y  over  [a,b]  to  be 

b  b 

(8.2)  (IR)Jy(x)dx  *  (R)Jy(x)dx, 

a  a 

provided  this  limit  exists  (in  the  extended  real  number  system;  infinite  values 
will  be  accepted  here  for  improper  integrals) .  Similarly,  if  y  is  Riemann  in- 
tegrable  over  [a,b]  for  b  >  a  finite,  then 


«  . .  b 

(8.3)  (IR)/y(x)dx  -  ££  (R)/y(x)dx 

a  a 

by  definition,  again  if  the  indicated  limit  exists. 

The  definition  of  interval  integrals  given  in  §6  yields  values  of  certain 
improper  Riemann  integrals  if  the  integrand  y  is  interpreted  to  be  the  degenerate 
interval  function  [y  ,  y] ,  for  example, 


(a) 

(8.4)  (b) 

(c) 

In  the  above,  the  value 


f1  -1/3.,  r3  , 

/  x  '  dx  *[-,«>) , 
0  * 


/  x  *dx  *  (“,*], 
0 


/  (-e~X)dx  -  [-»  ,  -i] . 

0 

of  the  improper  Riemann  integral  appears  as  the  finite 
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endpoint  in  each  of  the  intervals  (8.4) (a)  and  (8.4)  (c).  Tne  indegenerate  in¬ 
terval  (8.4) (b)  indicates  correctly  that  the  corresponding  improper  Riemann  in¬ 
tegral  is  divergent. 

Definition  8.1.  An  interval  integral  (6.1)  io  said  to  be  infinite  if  its 
value  is  one  of  the  indegenerate  intervals  [-«  ,  -°o]  or  [»,»],  indeterminant  if 
it  is  equal  to  R  *  [-<*>  ,  =°]  ,  or  Improper  if  its  value  is  a  semi-infinite  inter¬ 
val  [a  ,  »]  or  [-“  ,  b] ;  otherwise,  it  is  said  to  be  finite . 

The  relationship  between  improper  interval  and  Riemann  integrals  will  now 
be  considered  for  the  cases  (8.2)  and  (8.3),  as  illustrated  by  (8.4) (a)  and 

(8.4)  (c),  respectively. 

Suppose  that  y(x)  is  unbounded  above  at  x  =  a.  Thus,  every  Darboux  sum 
(4.7)  will  contain  a  term  of  the  form  (after  elimination  of  nondistinct  partition 
points,  if  necessary), 

(8.5)  w(X1)*Vy1=  [wtt^  «c  ,  *]  , 


where  X1  =  [a  ,  x^  ,  and 


(8.6) 


°i  "  xett.ty(x,K 


The  interval  integral  of  y  will  hence  either  be  improper  or  infinite.  The  fol¬ 
lowing  theorem  is  illustrated  by  (8.4) (a). 

Theorem  8.1.  Suppose  that  y  is  Riemann  integrable  over  [a  ,  b]  for  a  <  a  £  b, 
and  the  indefinite  interval  integral  I  y(ct)  satisfies 

cl 


(8.7) 


lim 

a -I 


“  Vl8)  =  aia  ^y(x)dx  =  10 


then,  the  improper  Riemann  integral  (8.2)  of  y  over  [a,b]  exists,  and 


(8.8) 


b  b 

/ y (x) dx  =  [ (IR) fy (x)dx  ,  «] . 
a  a 


Proof :  One  has 

bob 

(8.9)  fy (x) dx  =  / y(x)dx  +  / y(x)dx, 

a  a  a 


by  Theorem  6.3,  and,  by  Remark  6.3, 


(8.10) 


b  b  b  .  b 

/ y(x)dx  ■  [(R) Jy(x)dx  ,  (R)/y(x)dx]  «  (R)/y(x)dx, 
a  a  a  a 

as  degenerate  intervals  may  be  identified  with  the  corresponding  real  numbers. 

Taking  the  limit  as  a+a  of  both  sides  of  (8.9)  gives  (8.8).  QED. 

In  the  case  of  integration  over  an  infinite  interval,  say  [a  ,  »] ,  suppose, 

for  example,  that  y  is  negative,  but  that  y(x)t0  as  x  +  »,  as  in  (8.4) (c).  Then, 

each  Darboux  sum  (4.7)  will  correspond  to  a  partition  A  with  x  .  finite,  x  = 

n  n-1  n 


+®,  and  as 

(8.11) 


Vyn  =  [c  ,  o] , 


where  c  =  c  (y(x)}  <  0,  w(X  )  =  w([x  ,<»])=  ®,  then  each  will  contain  a  term 
xtx  n  n-1 

n 

equal  to 

(8.12)  w(X  )«W  =  [-»  ,  0], 

n  n 

by  rules  (2.1) (iii)  and  (2.1) (ix) .  The  situation  illustrated  by  the  example 

(8.4) (c)  is  a  case  of  the  following  result. 

Theorem  8.2.  Suppose  that  y  is  Riemann  integrable  over  the  finite  interval 

00 

[a,b]  for  each  b  >  a,  and  the  indefinite  interval  integral  I  y(b)  satisfies 

(8.13)  I°y(b)  =  [-»  ,  0); 

then,  the  improper  Riemann  integral  (8.3)  of  y  over  [a  ,  ®]  exists,  and 

(8.14)  /y(x)dx  =  I-®  ,  (IR)  Jy(x)dx] . 

a  a 

Proofs  This  follows  exactly  in  the  same  way  as  Theorem  8.1  by  writing 


(8.15) 

and  noting  that 

(8.16) 


/ y(x)dx  *  / y(x)dx  +  /y(x)dx, 
a  a  b 


b  b 

/y(x)dx  »  (R)/y(x)dx 

a  a 


as  a  degenerate  interval.  QED. 

Other  cases  of  improper  interval  and  Riemann  integrals  may  be  treated  in 
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a  similar  fashion. 

9.  Computational  implications  of  the  theory.  One  purpose  of  the  theory  of 
integration  of  interval  functions  developed  above  is  to  provide  a  theoretical 
framework  for  the  investigation  of  the  numerical  solution  of  linear  and  nonlinear 
integral  equations  such  as 

b 

(9.1)  u(x)  =  (R) Jg (x,t,u(x) ,u(t) )dt 

a 

by  interval  methods.  One  approach  along  these  lines  is  to  reformulate  (9.1)  as 


an  interval  equation. 


U  -  T(U) , 


for  an  interval  function  U  which  contains  the  desired  solution  u  of  the  integral 
equation  (9.1).  Under  certain  conditions,  the  operator  T  will  be  a  contraction 
mapping  ([1],  [2]),  and  the  iteration  process 


(9.3) 


Un+1  *  T(U^),  n  *  0,1,2,... 


will  converge  to  give  a  solution  of  (9.2) .  To  implement  this  for  the  integral 

equation  (9.1),  one  forms  the  interval  functions  G  "[g_,g],n»  0,1,2,..., 

n  *n  n 


where 


(9.4) 


and  then  (9.3)  becomes 


(9.5) 


c^x.t)  -  inf{g(x,t,Un(x)  ,un(t))  }, 


gn(x,t)  ■  sup{g(x,t,Un<x)  ,Un(t))  }r 


U  A.(x)  -  Jg  (x,t,U  (x),U  (t) )dt 

Ti+ 1  '  11  n  ii 


in  terms  of  interval  integration,  of  course,  if  ^(x, t)  and  gn(x,t)  are  Riemann 
integrable  in  t,  then  the  endpoint  functions  Hn+^,  ur+^  of  are  obtained  by 

Riemann  integration.  From  a  numerical  standpoint,  in  this  case  approximations 
u*+1  <_  Bn+^*  u**i  _>  ur+^  may  be  obtained  to  prescribed  accuracy  by  any  one  of  a 
number  of  methods,  including  the  use  of  Darboux  sums  as  defined  in  14  [7],  with 
higher-order  accuracy  being  obtainable  from  integration  of  Taylor  polynomial  ap¬ 
proximations  to  the  endpoint  functions,  or  by  other  rules  of  numerical  integration 


( [4] ,  [5],  [6],  [9])/  provided,  of  course,  that  the  endpoint  functions  are  smooth 
enough. 


A  particularly  simple  case  occurs  if  g  is  monotone  in  the  sense  that 

g  (x,t)  =  g (x,t,u  (x),u  (t) ) , 

— n  -n 


(9.6) 


g  (x,t)  =  g (x,t,u  (x),u  (t) ) , 
n  n  n 

that  is,  the  endpoint  functions  of  transform  into  the  endpoint  functions  of  G^, 
and  furthermore,  g  transforms  Riemann  integrable  functions  into  Riemann  integra- 
ble  functions.  Here,  the  iteration  (9.5)  can  be  carried  out  using  only  the  end¬ 
point  functions  if  one  starts  with  an  interval  =  [u^  ,  uQ]  which  has  Riemann 
integrable  endpoint  functions.  An  example  of  this  approach  to  the  solution  of  a 
nonlinear  integral  equation  was  given  by  Rail  [7] ,  in  which  step-functions  were 
used  as  endpoint  functions  (and  T  was  approximated  by  a  numerical  operator  S  such 
that  T  c  s)  .  In  many  cases,  continuous  solutions  u  are  sought  for  integral  equa¬ 
tions  (9.1),  which  gives  rise  to  the  following  concept. 

Definition  9.1.  The  continuous  core  C  (U)  of  an  interval  function  U  on  [a,b] 
is  defined  to  be  the  set  of  continuous  functions  y  contained  in  U,  that  is 

(9.7)  CC(U)  =  {y  |  y  €  UHC[a,b] } . 

Evidently,  C  (U)  c  C  (U) ,  the  Riemann  core  of  u  defined  earlier. 

C  K 

If  g  is  a  continuous  function  of  its  arguments,  and  the  interval  operator 
T  is  such  that  the  continuous  function  v  defined  by 


b 

(9.8)  v(x)  *  (R) Jg(x,t,u(x) ,u(t) )dt 

a 

belongs  to  T(U)  for  u  €  CC(U) ,  then  it  follows  that  each  continuous  solution  u 
of  (9.1)  will  belong  to  C  (U)  if  U  is  a  fixed  interval  of  T,  that  is,  if  (9.2) 

v 

holds.  Thus,  it  is  tempting  to  try  to  compute  the  sequence  (9.3)  using  only 
Cc(Un),  where  UQ  is  taken  to  have  continuous  endpoint  functions.  However,  in 
general,  the  functions  c^(x,t)  and  g^(x,t)  obtained  from  (9.4)  will  be  only  semi- 
continuous  if  UQ  is  replaced  by  Cc(Uq),  and  these  so-called  L-  and  U-functions 
may  not  even  be  Riemann  integrable  13) .  The  theory  of  interval  integration 
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developed  in  this  paper  resolves  this  difficulty  by  allowing  computation  with 
the  interval  functions  directly,  regardless  of  the  character  of  their  endpoint 
functions. 

Remark  9.1.  If  u  €  C  (U  )  is  a  solution  of  (9.1),  then  for  the  sequence 

v  U 

(9.3)  constructed  by  the  operations  (9.4)  and  the  interval  integration  (9.5),  it 
follows  from  the  condition  (9.8)  for  continuous  g  that 

(9.9)  u  €  C  (U  ) ,  n  =  0,1,2,...; 

C  II 

furthermore ,  for 

00 

(9.10)  u  -  n  u  , 

i  n 

n=l 

one  has  u  €  Cc  (U) . 

Remark  9.2.  In  the  favorable  case  that  U  ,  c  u  ,  n  -  0,1,2, _ _  and 

n+1  n  ' 

u.u)  £  w»-o. 

one  has  that  U  =  [u  ,  ul  ■  u  defined  by  (9.10)  satisfies  the  integral  equation 

(9.1),  as  a  degenerate  interval  integral  of  a  degenerate  interval  function  is 

necessarily  a  Riemann  integral;  furthermore,  one  has  error  bounds  of  the  form 

(9.12)  u  (x)  <  u(x)  <  u(x),  a  <  x  <  b, 

— n  —  —  _  — 

for  n  *  0,1,2, ... . 

Further  applications  of  interval  integration  to  the  solution  of  integral 
equations  will  be  investigated  in  subsequent  papers. 
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function  and  Jy(x)dx  is  the  upper  Darboux  integral  of  the  upper  endpoint 

a 

function  y.  As  these  Darboux  integrals  always  exist  in  the  extended  real  number 
system,  it  follows  that  all_interval  functions  are  integrable,  no  matter  how  nasty 
the  endpoint  functions  y  are.  The  interval  integral  defined  in  this  way  in¬ 
cludes  the  interval  integral  of  R.  E.  Moore  as  the  special  case  that  y  are 
continuous,  and  hence  Riemann  integrable. 

In  addition  to  a  construction  of  the  interval  integral  in  a  form  suitable 
for  numerical  approximation,  some  of  its  basic  properties  and  other  implications 
and  applications  of  its  definition  are  presented.  The  theory  of  interval  inte¬ 
gration  given  here  supplies  a  previously  lacking  mathematical  foundation  for  the 
numerical  solution  of  integral  equations  by  interval  methods. 


